FREIB URG-EHEP-96-01 



1 



JET PRODUCTION IN PHOTON-PHOTON INTERACTIONS * 

STEFAN SOLDNER-REMBOLD 
Albert- Ludwigs-Universitdt Freiburg, Fakultdt fur Physik, D-79104 Freiburg im Breisgau, Germany 

On behalf of the OPAL collaboration 

The inclusive one- and two-jet cross-sections are measured in collisions of quasi-real photons at e+e - centre-of- 
mass energies y/s ea = 130 and 136 GeV using the OPAL detector at LEP. Jets are reconstructed with a cone 
jet finding algorithm. The jet cross-sections are compared to next-to-leading order (NLO) perturbative QCD 
calculations. Transverse energy flows in jets are studied separately for direct and resolved two-photon events. 



1 Introduction 

We present measurements of jet production in 77 
interactions at y/s ec of 130 and 136 GeV. The jet 
cross-sections are compared to predictions of per- 
turbative QCD which use different parametrisa- 
tions of the photon structure function. 

In the Quark Parton Model (QPM) jets are 
produced by the interaction of bare photons, 77 — > 
qq. This is called the direct process. The largest 
part of the total cross-section, however, is mod- 
elled by interactions where the photon fluctuates 
into a hadronic state. The processes are called 
single-resolved if one photon couples directly to a 
parton in the other photon and double-resolved if 
partons from both photons interact. 

2 Data analysis and 
Monte Carlo simulation 

Two-photon events are selected by requiring that 
the sum of all energy deposits in the calorime- 
ters is less than 50 GeV. The hadronic invariant 
mass, WecaL) measured in the electromagnetic 
calorimeter has to be greater than 3 GeV; the 
missing transverse energy of the event has to be 
less than 5 GeV; at least 5 charged tracks must 
have been found and no track may have a momen- 
tum greater than 15 GeV/c. 

To remove events with scattered electrons in 
the forward detectors, the energy measured in 
these detectors has to be less than 40 GeV (20 
GeV) depending on the angular range, corre- 
sponding to an effective maximum photon virtual- 
ly, Qmaxj of 0.8 GeV 2 ("anti-tagging condition"). 

For an integrated luminosity of 4.9 pb _1 , 7808 
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events remain after all cuts. The average vis- 
ible hadronic invariant mass W vls is about 18 
GeV. The total background from other processes 
is about 1 %. 
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Figure 1: The number of two-jet events as a function of 
x+ compared to PHOJET and PYTHIA. Statistical errors 
only are shown. The hatched histogram is the direct con- 
tribution to the PYTHIA events. 

The leading, order (LO) QCD generators 
PYTHIA 5.721 1 and PHOJET 1.05 are used 
for the 77 event simulation. The fragmentation 
is handled by JETSET 7.4080. 

3 Jet reconstruction and backgrounds 

In the cone jet finding algorithm a jet is defined 
as a set of particles, i. e. partons or hadrons gen- 
erated in a Monte Carlo program, or tracks and 
calorimeter clusters, whose momenta lie within a 
cone of size R, such that the axis of the cone co- 
incides with the momentum sum of the particles 
contained. The cone size i?, which is chosen to be 
R = 1, is defined as R = ^(Ary) 2 + (A</>) 2 , with 
X] = — lntan(0/2) being the pseudorapidity and <fi 
the azimuthal angle in the laboratory frame. Arj 
and Acj) are the differences between the cone axis 
and the particle direction. The total transverse 
energy of the jet (defined with respect to the 



beam axis) is the scalar sum of the transverse en- 
ergies of its components. The pseudorapidity 
rf ct and the azimuthal angle <^ et are defined as 
the sum over the rji and <f>i of the jet components 
weighted by their transverse energies ■ 

All results are given for jets with E^ > 3 GeV 
and \if et \ < 1. 

4 Energy flow in jet events 

In 7p scattering, direct and resolved two-jet 
events, which correspond to the single- and 
double-resolved events in 77 scattering if the pro- 
ton is substituted in the place of a VMD-like pho- 
ton, have been identified by measuring the fraction 
x~f of the nhoton energy participating in the hard 
scattering!! In 77 scattering two photons interact, 
therefore a pair of variables is defined □ 



hadr 



(i) 



where p z is the momentum component along the 
z axis of the detector and E is the energy. Ide- 
ally, the total energy of the event is contained 
in the two jets for direct events without remnant 



jets, i. e. 



1, whereas for resolved events at 



least one of the x^ values is expected to be much 
smaller than one. The variables x^ are measured 
using all tracks and calorimeter clusters that were 
used in the jet finding algorithm. In addition, the 
energies in the forward detectors, which are not 
used for jet finding and for the energy flow distri- 
butions, are added in the denominator of Eq. [I]. 
This improves the separation of direct and re- 
solved events. 

Figure [l] shows the number of two-jet events 
as a function of arf", together with the PYTHIA 
and PHOJET samples, after the detector simula- 
tion. The x~ distribution, which is not shown, is 
consistent with the x+ distribution within the sta- 
tistical errors. The peak expected for direct events 
at x^T = 1 is smeared out due to higher order 
QCD, hadronisation and detector resolution ef- 
fects. Only statistical errors are shown, since sys- 
tematic effects largely cancel in the ratio (Eq. [I]) . 

The direct events, which are shown sepa- 
rately for the PYTHIA sample, mainly contribute 



O 16 

14 

f 12 

g 8 

S 6 

?C 4 
H 

2 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

: OPAL <a) E 




, 1 1 1 | 1 1 1 1 

: OPAL 


1 1 1 1 | 1 1 1 
, — ♦— Data _i 


; min(x* x^)<0.8 




; min(x*,x^)>0.8 J 


J PHOJET : 

| PYTHIA - 


+ : 






\ i 












12-2-1012 

At)' At)' 



- 

lOPAL 

: min(x* xl.)>0.8 


1 ' * 1 ■ t 

■ ] 


{ 

K..J 





Figure 2: Transverse energy flow measured relative to the 
direction of the jet in two-jet events. Jets from events 
with m'm(x^,Xy) < 0.8 (a,c) and min(xj, x~^ ) > 0.8 (b,d) 
are shown separately. The energy flow is integrated over 
|A</>| < it/2 for the A?/ projection (a,b) and over |A?/| < 1 
for the A<f> projection (c,d). Statistical errors only are 
shown. 



in the region x+ > 0.8. Experimentally, sam- 
ples with large and small direct contributions are 
separated by requiring mia(x^ ,x^) > 0.8 and 
, x~ N 



min(:r+,:c~) < 0.8, respectively. About 95 % of 
all events in the region min(a;+,a;~) > 0.8 origi- 
nate from direct interactions in PYTHIA. 

Figure || shows the transverse energy flows 

1 dE T , 1 dE T 



A jct d(Ar/) 



and 



ATjet d(A0) 



(2) 



with respect to the jet direction for two-jet events 
with xnhx{x^,x~) < 0.8 (low x 7 ) and with 
min(x+,a;~) > 0.8 (high x 7 ), separately. The 
shape of the selected jet is seen clearly in both Ai]' 
and A(f>. No detector correction has been applied. 
The same tracks and clusters which are used for 
jet finding are included, with 



At/ = k(r)- rf ct ) and Acj) 



(3) 



The rapidity difference is multiplied by a factor 
k = ±1. The factor k is chosen event-by-event to 
be k = +1 for events with 2;+ > x~ and k = — 1 
for events with x^ < XZ. Due to the choice of k, 
there is always more energy flow at Ar/ < 0. 
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In Fig. Q, a pedestal is observed which could 
indicate the existence of a photon remnant. As 
expected, this enhancement in the region around 
|A0| Ri 7r/2 and at |Ar/| > 1 is more pronounced 
for low x 1 events. The jets in high x^events are 
much more back-to-back in A(j) (Fig. ||d) than in 
low x 1 events (Fig. ||c). The pedestal in the A0 re- 
gion between the two jets at | A</>| i=s 7r/2 is not ob- 
served in the high x 7 events. Jets in high cc 7 events 
are also observed, on average, to have more aver- 
age transverse energy and to be more collimated. 
This is as expected for direct events, where all the 
available energy is used in the hard subsystem. 
Both Monte Carlo models describe the transverse 
energy flow reasonably well, except for an under- 
estimate in the central region around the jet axis. 
Nevertheless the overall modelling is sufficiently 
good to justify using the Monte Carlo models for 
unfolding the detector resolution effects in the jet 
cross-section measurements. 

5 Inclusive jet cross-sections 

The one- and two-jet cross-sections shown in Fig. || 
are corrected to the hadron level using the method 
of regularised unfolding. The main systematic un- 
certainty on the jet cross-sections originates from 
varying the energy scale of the EC AL in the Monte 
Carlo simulation by ±5 % and from the unfolding 
procedure. The model dependence of the unfold- 
ing is taken into account by adding to the system- 
atic error the difference between the results ob- 
tained with PYTHIA, which are taken to be the 
central values, and PHOJET. 

The distribution is compared to a NLO 
perturbative QCD calculation of the inclusive one- 
jet cross-section by Kleinwort and Kramer a who 
use the NLO GRV patametrisation of the pho- 
ton structure function Their calculation was 
repeated for the kinematic conditions of this anal- 
ysis. All scales are chosen to be equal to E 1 ^ . 
The strong coupling a s is calculated from the two- 
loop formula with A^l = 130 MeV. The ratio 

MS 

of the NLO to the LO inclusive one-jet cross- 
section decreases from about 1.19 to 1.03 between 
E^ = 3 GeV and E^ = 16 GeV. 

The direct, single- and double- resolved parts 
of the one-jet cross-section and their sum are 
shown separately. The agreement between data 
and the calculation is good. The resolved cross- 



sections dominate in the region £Jf < 5 GeV, 
whereas, at high E^ the direct cross-section is 
largest. The E 3 ^ distribution falls less steeply 
than expected in 7p and pp interactions, because 
the fraction of hard interactions rises from pp to 
7P to 77 interactions, due to the additional direct 
photon interactions. It should be noted that the 
NLO QCD calculation gives the jet cross-section 
for massless partons, whereas the jet cross-sections 
are measured for hadrons. 

The inclusive two-jet cross-section is measured 
for events with at least 2 jets using the two jets 
with the highest E^ in an event. The E^ 1 distri- 
bution is also compared to the NLO calculationsEI, 
which have been extended to include the full NLO 
calculation for the double-resolved two-jet cross- 
section. The increase of the NLO compared to the 
LO inclusive two-jet cross-section is less than 10 % 
in the calculation. 

The inclusive one- and two-jet cross-sections 
as a function of |??i et | are shown in Fig. |J. Within 
the uncertainties of the measurement, the distri- 
butions are nearly independent of \if et \ in agree- 
ment with the expectations of the Monte Carlo 
models. 

The total cross-sections, which are dominated 
by the low Ej? events, depend on the photon 
structure function. In Fig. ^ the total jet cross- 
sections predicted by the two Monte Carlo mod- 
els differ significantly even if the same photon 
structure function (here SaS-lD) is used. This 
model dependence reduces the sensitivity to the 
parametrisation of the photon structure function. 
Different parametrisations were used as input to 
the PHOJET simulation. The GRV-LO and SaS- 
lD parametrisations describe the data equally 
well, but the LAC1 parametrisation overesti- 
mates the total jet cross-section by about a factor 
of two. 

6 Conclusions 

We have measured jet production in 77 interac- 
tions with the OPAL detector at v^oc °f 1^0 an d 
136 GeV, using a cone algorithm with R = 1, 
E^ > 3 GeV and |ryJ et | < 1. 

Two-jet events originating from direct and 
resolved photon interactions were separated ex- 
perimentally using x^ . Jets in events with 
min(a; 7 ", x~) > 0.8 are expected to be produced 
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Figure 3: The inclusive one-jet (upper plot) and two-jet 
(lower plot) cross-sections as a function of E}£ .compared 
to the NLO calculation by Kleinwort and KramerQ. The di- 
rect, single-resolved and double-resolved cross-sections and 
the sum are shown separately. 
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Figure 4: The inclusive one-jet (upper plot) and two-jet 
(lower plot) cross-sections as a function of l^ '! compared 
to the LO QCD calculations of PYTHIA and PHOJET. 



mainly from direct interactions. These jets are 
observed to have, on average, more average trans- 
verse energy and to be more collimated than 
jets in resolved events with min(2;+, £~) < 0.8. 
In resolved events a pedestal is observed in the 
transverse energy flows which may be related to 
the photon remnant. The Monte Carlo models 
PYTHIA and PHOJET describe the transverse 
energy flow distributions reasonably well. 

The inclusive one-jet and two-jet production 
cross-sections were measured as a function of E^ 
and \if ct \. The measurement extends the Eif? 
range of previous measurements u up to Ej, = 
16 GeV. The E 1 * 1 dependent one- and two-jet 
cross-sections are .in good agreement with NLO 
QCD calculations 

Within the uncertainties of the measurements, 
the jet cross-sections are nearly independent of 
\rf et \ for \if et \ < 1. The total jet cross-section 
is dominated by the resolved cross-section. Given 
the model dependence of the jet cross-sections in 
PYTHIA and PHOJET, the GRV-LO and SaS-lD 
parametrisations describe the data equally well. 
The LAC1 parametrisation overestimates the to- 
tal jet cross-section by about a factor of two. 
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